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I.  INTRODUCTION 


The  Euler  Tunnel  Analysis  (ETA)  suite  has  been  demonstrated  to  be  a  valuable  tool  in  the 
aerodynamic  design  of  missiles  at  the  U.S.  Army  Aviation  and  Missile  Research,  Development 
and  Engineering  Center  (AMRDEC)  [1].  It  provides  a  high-fidelity,  complementary  collaborator 
to  quick,  preliminary  design  methods  such  as  Missile  DATCOM  [2],  AP05  [3],  and  Missile3  [4], 
And  it  offers  a  straightforward,  easy-to-use  alternative  to  the  fragmented  approach  of  employing 
multiple  intermediate-  to  high-fidelity  techniques  (like  the  S/HABP  [5]  piecemeal  method;  the 
CMARC  [6]  and  PANAIR  [7]  panel  codes;  the  ZEUS  [8]  marching  technique;  parabolized 
Navier-Stokes  schemes;  and  full-field,  elliptical  approaches,  that  is,  potential  flow,  Euler,  and 
Navier-Stokes  Computational  Fluid  Dynamics  (CFD)  solvers)  to  construct  the  flight 
characteristics  of  a  missile. 

ETA  does  this  by  assembling  the  components  needed  to  conduct  an  Euler  flowfield  analysis 
(including  an  automated  grid  generator  and  a  robust  solver)  within  a  single  framework.  This 
substantially  reduces  the  time  and  specialized  knowledge  required  to  generate  grids  and 
setup/operate  the  solver,  resources  that  are  often  unavailable  to  the  aerodynamic  designer.  In 
addition,  the  ETA  solver  incorporates  the  vorticity  confinement  technique  [9]  to  conserve  field 
vortices  and  mimic  surface  turbulence.  With  these  attributes,  ETA  is  capable  of  greater 
productivity,  physical  fidelity,  and  accuracy  than  is  typical  of  most  Euler  approaches. 

II.  METHODOLOGY 

A.  Description 

ETA  is  a  suite  of  Government-owned,  productivity-oriented,  CFD  software  developed 
to  facilitate  aerodynamic  design  and  analysis.  As  mentioned  above,  it  contains  an  automated  grid 
generator  and  a  robust  Euler  solver.  However,  it  also  contains  a  geometry  generation  code, 
named  CFDGEN,  to  construct  configurations  from  a  library  of  pre-existing  models  and  model 
parts.  CFDGEN  does  not  require  Computer  Aided  Design  /  Computer  Aided  Manufacturing 
(CAD/CAM)  expertise  so  the  aerodynamic  designer  can  use  it  to  construct  body  geometries 
without  specialized  training.  Although  straightforward  to  use,  CFDGEN  is  not  visually 
interactive  and  thereby  cumbersome  to  apply.  To  circumvent  this  inconvenience,  the  AMRDEC 
often  utilizes  some  of  the  less  expensive  commercial  CAD/CAM  packages  to  construct,  convert, 
and  repair  configurations  of  interest. 

With  regards  to  automated  grid  generation,  ETA  exercises  the  cubes  component  of  the 
NASA/Ames  Cart3D  [10]  methodology  developed  by  Aftosmis,  Melton,  and  Berger.  This 
technique  produces  an  unstructured,  Cartesian  field  grid  around  the  body  with  relatively  few  user 
inputs.  The  time  required  to  do  so  is  typically  less  than  an  hour  for  most  of  the  AMRDEC ’s 
cases — an  attribute  that  saves  numerous  man-hours  when  setting  up  computations.  In  addition, 
ETA  utilizes  the  Cart3D  geometry  import  and  conversion  programs.  Hence,  it  interfaces  with 
several  geometry  formats.  However,  the  NASA  Langley  Wireframe  Geometry  Standard 
(LaWGS)  format  (produced  by  CFDGEN)  and  the  stereolithography  (STL)  solid  model  format 
(produced  by  many  CAD/CAM  packages)  are  the  ones  most  often  used  by  the  AMRDEC. 


The  flow  solver  component  of  ETA  is  the  NASA/Ames  Topology  Independent 
Gridding  Euler  Refinement  (TIGER)  [11]  code  originally  developed  by  Melton  and  modified  by 
Robinson  [12,13].  It  is  an  explicit,  unstructured,  hexahedral,  finite-volume,  Euler  CFD  method 
that  employs  Jameson’s  Runge-Kutta  scheme  [14]  for  time  integration.  The  enhancements  made 
by  Robinson  extend  the  code  by:  (1)  permitting  use  of  an  algebraic  enthalpy  equation  in  lieu  of 
the  differential  energy  equation,  (2)  generalizing  the  Runge-Kutta  integration  scheme  from  its 
original  four-stages  to  m-stages  specified  by  the  user  (usually  two),  (3)  enhancing  robustness  near 
corners  and  other  high- gradient  regions  at  high  Mach  numbers,  and  (4)  adding  the  vorticity 
confinement  technique  [9]  to  conserve  vorticity  in  the  field  and  over  solid  bodies. 

A  Graphical  User  Interface  (GUI)  and  several  software  scripts  enable:  (1)  editing  of 
body  geometries  and  the  deflection  of  control  surfaces,  (2)  generation  of  the  solution  grid, 

(3)  establishment  of  a  directory  structure  and  flow  solver  input  files  (based  on  Mach  number,  roll 
angle,  and  angle  of  attack),  (4)  creation  of  the  required  directory  structures  and  input  files, 

(5)  submission  of  cases  to  the  computational  platform,  (6)  checking  the  status  of  each  case, 

(7)  retrieval  of  the  results,  and  (8)  calculation  of  aerodynamic  coefficients  from  each  flowfield 
solution  for  each  point  of  interest. 

B.  Initial  Application 

ETA  was  initially  applied  to  an  eight-finned  missile  configuration  (Fig.  1)  that  uses  a 
bent  nose  to  effect  aerodynamic  control.  The  baseline  geometry  was  constructed  with  CFDGEN 
using  wind  tunnel  model  blueprints  and  it  incorporated  three-dimensional  fin  thickness  as  well  as 
edge  radius,  thickness  taper,  and  breakline  sweep  angles.  The  nose  angle  was  then  varied  to 
create  individual  representations  for  each  deflection  angle. 

For  each  bent  nose  geometry,  the  input  to  cubes  was  set  to  create  a  field  grid  with  the 
upstream  boundary  approximately  0. 19  body  lengths  away  from  the  missile.  The  downstream 
boundary  was  1  length  away,  and  the  side  boundaries  were  located  from  1.5  to  1.9  lengths  away 
(Fig.  2).  The  flow  conditions  were  set  to  freestream  values  on  each  boundary,  while  the 
minimum  grid  cell  dimensions  were  set  to  approximately  one  one-thousandth  (1/1000)  of  the 
body  length.  Six  levels  of  refinement  were  also  specified  which  resulted  in  grids  populated  with 
approximately  one  million  (1,000,000)  Cartesian  cells.  Computations  were  then  performed  in 
strictly  Euler  mode  (without  vorticity  confinement)  for  a  Mach  number  of  3.0;  roll  angles  of  0, 

45,  and  90  degrees;  and  angles  of  attack  from  -10  to  +10  degrees. 
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Figure  1.  Bent-Nose  Missile  -  Delta  8  Degree  Nose 


Figure  2.  Typical  Grid  for  Bent-Nose  Missile  Bodies 
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C.  Surface  Confinement  Coefficient 


The  implementation  of  vorticity  confinement  in  ETA  is  described  in  Reference  13 
which  states,  “The  essential  feature  of  compressible  vorticity  confinement  involves  the  addition 
of  an  anti-diffusive  term  to  the  momentum  equations,  proportional  to  the  local  vorticity 
magnitude  and  in  a  direction  opposed  to  the  numerical  diffusion  of  vorticity  ...”  This  term  takes 
different  forms  for  field  and  surface  corrections,  but  each  approach  requires  a  user-specified 
input  constant  to  establish  the  degree  of  correction  that  will  be  applied.  Experience  has 
demonstrated  that  the  default  field  confinement  coefficient  in  ETA  (established  by  Robinson) 
appears  to  be  adequate  for  the  most  part.  However,  the  surface  confinement  coefficient  is  more 
problematic  since  the  corrections  it  governs  are  directly  related  to  viscous  flow  phenomena 
(usually  turbulent)  at  and  near  the  wall  for  both  attached  and  separated  flows.  Since  there  are 
currently  no  general,  physically-based  guidelines  for  determining  appropriate  values  of  this 
constant,  an  attempt  was  made  to  rectify  this  situation  for  the  problem  at  hand. 

If  the  Euler  equations  with  surface  confinement  terms  are  equated  to  the  Favre- 
averaged  form  of  the  Navier-Stokes  equations,  a  relationship  can  be  established  between  the 
confinement  parameter  and  the  viscous  terms.  Then  the  approach  used  in  developing  the 
compressible  Mentor  Shear  Stress  Transport  (SST)  turbulence  model  [17,  18]  is  followed.  In 
this  method,  the  turbulence  production  term  is  presumed  to  be  identical  for  both  compressible 
and  incompressible  forms  of  the  Navier-Stokes  equations.  Since  the  flow  adjacent  to  the  wall  is 
very  likely  to  be  incompressible  in  the  laminar  sublayer,  the  incompressible  form  of  the  equations 
is  used.  By  assuming  turbulence  production  to  be  homogeneous  and  turbulent  viscosity  to  be 
spatially  constant  for  the  local  flow,  it  is  possible  to  express  the  production  in  terms  of  vorticity. 
Since  the  laminar  and  turbulent  viscous  terms  of  the  Navier-Stokes  equations  have  the  same  form 
(due  to  the  Boussinesq  assumption)  the  laminar  components  of  the  equations  can  likewise  be 
articulated  in  terms  of  vorticity. 

The  net  result  is  (utilizing  Robinson’s  notation  in  Reference  13,  Anderson  et.  al.’s 
notation  in  Reference  19,  and  Hinze’s  notation  in  Reference  20) 


£s  ph2 


(1) 


where  £s  is  the  surface  confinement  coefficient,  p  is  the  fluid  density,  h  is  the  representative 
dimension  of  the  computational  cell,  w  is  the  local  vorticity,  //  =  jUlairiaar  +  plurhuk,nl  (he.,  the  sum 
of  the  laminar  and  turbulent  eddy  viscosities),  xj  is  the  Cartesian  coordinate  in  the  j  direction, 
and  Uj  is  the  local  velocity  component  in  the  i  direction.  By  performing  an  order  of  magnitude 
analysis  on  the  various  terms  this  equation  simplifies  to 
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(2) 


which  relates  the  surface  confinement  coefficient  to  a  combined  laminar/turbulent  form  of 
Reynolds  number. 

From  Wilcox  [21]  the  “outer”  part  of  the  Cebeci-Smith  “mixing  length”  turbulence  model  is 
shown  to  be  given  by 


=0'0168/4,mi[W  Rer 


(3) 


(where  Re^.  is  the  Reynolds  number  based  on  displacement  thickness)  so  that  the 
coefficient  can  now  be  written  in  a  general  form  given  as 


(4) 


puh 


As  a  test,  this  formula  was  exercised  for  flow  over  a  turbulent  flat  plate  in 
incompressible  flow.  It  was  modified,  though,  by  replacing  the  displacement  thickness  with 
boundary  layer  thickness.  Although  this  was  a  rough  engineering  approximation,  comparisons 
with  Dadone’s  corresponding  case  [22]  revealed  that  the  modified  formula  produced  a  surface 
confinement  coefficient  quite  similar  to  those  used  by  Dadone  et.  al. 

For  the  problem  at  hand,  Reference  23  was  used  to  find  the  boundary  layer  thickness  at 
the  separation  point  of  a  cylinder  in  cross  flow  and  relate  it  to  the  cylinder  radius.  Although  it 
was  only  possible  to  find  the  thickness  of  a  laminar  boundary  layer,  the  ratio  of  measured 
turbulent  and  laminar  boundary  thicknesses  in  Reference  24  was  used  to  approximate  a  turbulent 
layer.  In  addition,  since  displacement  thickness  is  a  more  difficult  parameter  to  obtain  than  layer 
thickness,  the  displacement  thickness  was  replaced  with  boundary  layer  thickness,  as  before.  The 
consequence  of  all  these  manipulations  was 


(5) 


which  casts  the  coefficient  in  terms  of  Reynolds  number,  cylinder  diameter,  and  cell  size. 

This  equation  was  applied  to  the  missile  airframe  but  produced  a  coefficient  that  was 
nearly  zero.  From  reviewing  the  literature,  it  was  judged  that  this  was  unworkably  small  and  that 
a  compressibility  correction  was  needed.  After  some  experimentation,  it  was  decided  to  recast 
Equation  5  as 


(6) 
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for  flow  scenarios  when  the  freestream  Mach  number,  MrM ,  is  greater  than  1 .  When  employed 

for  the  bent  nose  missile,  Equation  6  yielded  a  coefficient  within  4  percent  of  the  default  value. 
Since  the  calculated  coefficient  rounded  to  default  value,  the  default  setting  was  used  for  the 
surface  confinement  calculations  that  followed. 

D.  Application  with  Surface  Confinement 

To  test  the  veracity  of  Equation  6,  ETA  was  reapplied  to  the  missile  with  an  8-degree 
nose  deflection.  This  airframe  was  selected  because,  as  is  shown  in  the  subsequent  section,  it 
produced  the  greatest  disparity  between  computations  and  measurement.  For  these  calculations, 
the  default  grid  settings  were  used  so  the  downstream  boundary  was  nearly  two  body  lengths 
from  the  base,  and  the  other  bounds  were  about  one  length  away.  Again,  the  minimum  cell 
dimensions  were  set  to  roughly  one  one-thousandth  (1/1000)  of  the  body  length.  Nine  levels  of 
refinement  were  specified,  resulting  in  a  grid  composed  of  less  than  900,000  cells.  Computations 
were  then  performed  for  a  Mach  number  of  3.0;  roll  angles  of  0, 45,  and  90  degrees;  and  angles 
of  attack  ranging  from  -10  to  +10  degrees.  The  same  grid  was  used  for  all  these  calculations 
which  were  made  without  and  with  both  field  and  surface  confinement.  Although  there  are  no 
upstream  control  surfaces  to  generate  tip  vortices,  the  deflected  nosetip  could  create  a  vortex  that 
affects  the  body  flow,  especially  when  it  is  rolled  and  at  a  negative  angle  of  attack.  So  the  field 
mode  was  applied  as  appropriate.  Note  that  the  default  constants  were  utilized  for  both  modes  of 
confinement. 

III.  RESULTS  AND  DISCUSSION 

A.  Initial  Application 

Results  for  the  initial  application  of  ETA  to  the  bent  nose  configuration  are  presented 
in  Figures  3  through  6.  They  are  shown  in  the  wind  tunnel  reference  frame  and  compared  with 
measurements  documented  in  Reference  16.  The  predictions  were  performed  “blind”  without 
examining  any  of  the  wind  tunnel  data.  Yet  it  is  evident  from  Figures  3  and  4  that  the  ETA 
values  for  normal  force  and  pitching  moment  agree  well  with  the  data  for  each  of  the  bent  nose 
deflections  (each  at  a  roll  angle  of  0  degrees).  In  particular,  the  pitching  moment  agrees  very 
well  especially  since  it  is  referenced  to  the  center-of-gravity  (a  position  that  amplifies  prediction 
discrepancies).  It  should  also  be  observed,  though,  that  the  pitching  moment  predictions  begin  to 
deviate  from  the  measurements  when  the  sum  of  the  nose  deflection  and  the  angle  of  attack 
begins  to  exceed  10  degrees.  This  behavior  is  expected  since  viscous  separation  is  likely  to  occur 
above  moderate  angles  of  incidence. 
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Figure  3.  Normal  Force  Coefficient  Versus  Angle  of  Attack  for  All  Bent-Nose  Configurations 


Angle-of-Attack,  Deg 


Figure  4.  Pitching  Moment  Coefficient  Versus  Angle  of  Attack  for  All  Bent-Nose  Configurations 
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At  roll  attitudes  of  45  and  90  degrees,  the  normal  force  and  pitching  moment  curves 
were  nearly  coincident  with  the  measurements,  so  these  curves  are  not  shown  to  save  space. 
Rather,  the  more  demanding  comparison  of  yawing  moment  is  provided  in  Figures  5  and  6.  It 
can  be  seen  that  the  agreement  with  measurement  is  again  quite  good  up  to  about  4  degrees  angle 
of  attack.  Above  this,  as  mentioned  before,  viscous  separation  is  likely  to  have  occurred. 


Figure  5.  Yawing  Moment  Coefficient  Versus  Angle  of  Attack  for  All  Bent-Nose 

Configurations  at  45-Degree  Roll 


Figure  6.  Yawing  Moment  Coefficient  Versus  Angle  of  Attack  for  All  Bent-Nose  Configurations 

at  90-Degree  Roll 
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B.  Application  with  Surface  Confinement 

Products  of  the  reapplication  of  ETA  with  vorticity  confinement  to  the  8  degree  bent 
nose  airframe  are  presented  in  Figures  7  through  10.  In  Figure  7,  it  can  be  seen  that  vorticity 
confinement  improves  the  calculated  normal  force  curve  to  the  point  that  it  is  indistinguishable 
from  the  wind  tunnel  data.  And  in  Figure  8,  it  can  be  observed  that  between  -4  and  0  degrees  of 
alpha,  the  technique  does  not  agree  with  the  pitching  moment  data  as  well  as  the  strictly  Euler 
results.  However,  it  betters  those  results  by  capturing  the  trend  of  the  measurements  between  4 
and  10  degrees  where  flow  separation  occurs. 

The  rolled  body  outcomes  for  yawing  moment  are  exhibited  in  Figures  9  and  10.  Both 
of  these  plots  reveal  that  while  the  confinement  correction  produces  a  slight  bias  in  moment 
between  ±4  degrees,  it  also  begins  to  capture  the  correct  “roll  off’  in  moment  that  occurs  outside 
these  bounds.  Although  the  favorable  comparison  with  measurement  is  not  as  comprehensive  as 
desired,  it  is  clear  that  the  correct  kind  of  flow  physics  is  being  approximated.  In  short,  it  is 
evident  that  the  confinement  technique  is  able  to  account  for  flow  separation  on  the  smooth, 
cylindrical,  missile  body  whereas  the  strictly  Euler  approach  is  not.  This  is  a  significant 
development  in  the  application  of  Euler  methods.  In  addition  to  these  observations,  it  is  also 
evident  that  even  though  the  surface  confinement  constant  may  need  some  adjusting,  the  value 
provided  by  Equation  6  was  appropriate.  This  lends  credence  to  its  validity.  It  may  be  that  with 
some  further  modification,  this  formula  could  provide  a  constant  that  would  enable  ETA  to 
match  all  the  yawing  moment  data  shown  in  these  figures. 


Figure  7.  Normal  Force  Curve  for  Delta  8  Degree  Nose  -  With/Without  Confinement 
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Figure  8.  Pitching  Moment  Force  Curve  for  Delta  8  Degree  Nose  -with/Without  Confinement 
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Figure  9.  Yawing  Moment  Curve  for  Delta  8  Degree  Nose  at  45-Degree  Roll  - 

With/Without  Confinement 
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Figure  10.  Yawing  Moment  Curve  for  Delta  8  Degree  Nose  at  90-Degree  Roll  - 

With/Without  Confinement 


IV.  SUMMARY 

It  has  been  shown  that  the  vorticity  confinement  technique  enables  an  Euler  solver  to 
account  for  the  aerodynamic  effects  of  crossflow  separation.  This  was  demonstrated  in  particular 
for  a  missile  characterized  by  a  high  fineness  ratio  body  with  a  deflected  nose  at  roll  angles  of  0, 
45,  and  90  degrees.  In  addition,  a  relation  has  been  derived  that  connects  (in  general  terms)  the 
surface  confinement  coefficient  to  physical  parameters.  When  applied  to  the  circular  missile 
body  in  crossflow,  the  coefficient  produced  by  the  relation  was  found  to  yield  physically 
meaningful  results.  In  other  words,  the  ensuing  normal  force,  pitching  moment,  and  yawing 
moment  values  were  all  improved  significantly  over  those  produced  by  strictly  Euler 
computations.  And  this  enhancement  was  obtained  at  the  low  cost  of  modestly  increased 
computer  run  time  for  each  Euler  flow  calculation.  Although  the  improvement  was  not  as 
comprehensive  as  desired,  it  is  anticipated  that  further  investigation  will  yield  a  more  complete 
form  of  the  equation. 
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